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Рассмотрены математические аспекты решения пространственной динамической задачи контактирования 
цилиндрических деталей машин. Приведена модификация метода суперпозиции, позволяющая провести 
аналитическое исследование напряженно-деформированного состояния и дать практические рекомендации 
по оценке интенсивности локальной концентрации напряжений в нерегулярных зонах контакта, что предо- 
пределяет прочностные параметры контактной пары. В качестве численного эксперимента рассмотрены ме- 
ханические характеристики взаимодействия поршневых пальцев двигателей внутреннего сгорания с голов- 
кой шатуна и бобышками поршня. 
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Введение. Детали, воспринимающие контактные нагрузки, как правило, являются одними из 
наиболее ответственных, и от их прочности зависит работоспособность и надежность всего меха- 
низма в целом. Однако при рассмотрении систем контактирующих тел механическими явлениями 
в зоне контакта зачастую пренебрегают, проводя упрощение и схематизацию усилий, действую- 
щих на детали механизма. Такие упрощения могут приводить к неверным результатам, так как в 
большинстве реальных механизмов закон распределения контактных давлений оказывает сущест- 
венное влияние на напряженно-деформированное состояние (НДС) механизма или конструкции в 
целом. В связи с этим возникает необходимость исследования НДС деталей механизма в местах 
их контактного взаимодействия при произвольной площадке контакта, размеры которой зависят 
от приложенных нагрузок. 

Теоретические исследования в этой области связаны с решением краевых задач неста- 
ционарной динамики со смешанными граничными условиями. Одним из наиболее эффективных 
является аналитический подход, предложенный в [1] и развитый в работе [2]. В соответствии с 
этим подходом решение упругопластической задачи получают на основе комбинации решений 
упругой и жесткопластической задач. Такой подход базируется на принципе Герстнера, при кото- 
ром общее перемещение можно представить в виде суммы упругой и пластической составляю- 
щих, рассматриваемых независимо. 

Упругое поведение большинства конструкционных материалов играет важную роль в де- 
формационном процессе, особенно когда зона пластических деформаций мала. Решение контакт- 
ных задач при упругопластическом деформировании в основном проводится на основе численных 
схем. Разработка аналитических подходов возможна при введении ряда гипотез [1,2] или эмпи- 
рических коэффициентов [3]. В последние годы широкое распространение при решении упруго- 
пластических задач получили численные методы, среди которых наиболее универсальным явля- 
ется метод конечного элемента. 

Несмотря на большое количество работ, посвященных данной тематике, исследователи 
вынуждены повсеместно моделировать одно из контактирующих тел бесконечной или полубеско- 
нечной областью, что дает возможность применения при аналитическом решении методики инте- 
гральных преобразований [4], понижающей размерность задачи. Кроме того, существенным уп- 
рощением следует признать тот факт, что неограниченное тело моделируется областью с прямо- 
линейной границей (полупространство, слой). Такой подход, хотя и дает возможность математи- 
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ческого анализа исходной краевой задачи, однако, может быть уместен лишь в случае малой об- 
ласти контакта по сравнению с размерами контактирующих тел (точечный контакт). При анализе 
НДС деталей машин и особенно деталей двигателей внутреннего сгорания (ДВС) такая разновид- 
ность контакта встречается крайне редко. Поэтому проблему построения и решения математиче- 
ской модели динамической пространственной задачи контактирования упругих цилиндрических 
деталей конечных размеров следует признать актуальной. 

В настоящее время разработано два подхода к решению граничных задач динамической 
теории упругости для тел конечных размеров. Один из них, метод однородных решений [5,6], на- 
шел применение в плоской задаче теории упругости, в теории тонких и толстых плит. Решение 
задачи находится с помощью однородных решений, которые являются интегралами основных 
уравнений теории упругости и удовлетворяют нулевым граничным условиям на части поверхно- 
сти тела, совпадающей с одной из координатных поверхностей. Во втором подходе, развитом в 
[6-8], решение задачи представляется в виде суперпозиции нескольких последовательных част- 
ных решений, обладающих конкретными свойствами симметрии. При этом предполагается, что 
поверхность упругого тела образована частями координатных поверхностей разных семейств в 
декартовой, цилиндрической или сферической системах координат. Метод суперпозиции наибо- 
лее приспособлен для решения двумерных динамических задач теории упругости, что существен- 
но снижает области его применения. 

Целью данной работы является построение алгоритма метода суперпозиции, свободного 
от предположений раздельного рассмотрения волновых движений с различными типами симмет- 
рии относительно срединной плоскости цилиндра, распространения этого метода на трехмерные 
области с криволинейной границей и применения разработанной схемы для численно- 
аналитического расчета полей контактных напряжений в поршневом пальце (ПП) ДВС. 
Постановка задачи. Пусть ПП, моделируемый конечным полым цилиндром, занимает в цилинд- 
рической системе координат область: 

р: -Н<2<0\/Ь<В<а0<ф<2л. 

Предполагаем, что материал цилиндра характеризуется упругими параметрами: 

р (плотность), д, И (коэффициенты Ляме). Пусть вдоль границы области заданы следующие 





граничные условия в напряжениях( (2 -круговая частота вибронагружения): 
К=а: = фе, ть = фе, ть = И ф, ре: 
В=Ь: с =Л (ф,2)е, т, = (ф,2)е, „= (ф,2)е"; 
2=-Н:0 .=$ф (В, фе, ть. = (фе, т„=ф (К, фе"; 


2=0:0.=ф (В, фе", ть =$ (В, ф)е“, т. =ф (В, фе. 
Движение среды описывается уравнениями линейной теории упругости в смещениях, 
связь компонент векторов напряжений и деформаций определяется соотношениями закона Гука в 
цилиндрической системе координат [6]. 
Решение сформулированной краевой задачи строим с использованием принципа суперпо- 
зиции, используя тот факт, что данная область является пересечением подобла- 


стей 0, :а< К <Б и РО, :-Н <2< 0.Вектор смещения О(В,Ф,2) в области Др, 


(1) 


ищем в виде суммы векторов (У(К,ф,2)=и(К,ф,2) +и,(К,ф,2), определяющих решения 
вспомогательных задач для бесконечного цилиндрического слоя (2, - и) и слоя с плоскопа- 


раллельными границами ( Г), —и, ). 
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Начнем с решения вспомогательной задачи для бесконечного цилиндрического слоя. С 
этой целью рассмотрим установившиеся Гармонические колебания области 79) аль при 
следующих граничных условиях: 

. РВ + — |. + 14 М + 1“. 
К =а:0 „=А, (ф,2)е Ть. = Х› (ф,2)е То = Хз (ф,2)е > 


ыы —а — —а _ —0 (2) 
В=Ь: а „= АЛ, (ф,2)е ›Ть, = Х, (ф,2)е То = Хз (ф,2)е * 
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Функции напряжений Х 7, вообще говоря, неизвестны, но в рамках постановки вспомо- 


гательной задачи предполагаем их заданными и принадлежащими к классу суммируемых вдоль 
всей границы рассматриваемой подобласти. На бесконечности задаются условия излучения энер- 
гии колебаний. Эти предположения являются естественными и необходимы для корректного ис- 
пользования при решении задачи методики интегральных преобразований Фурье [4]. В результа- 
те достаточно сложных и громоздких выкладок получено следующее представление для компо- 


нент вектора перемещений: и(К,ф,2) = м = Инь И = Им» Из: = и, 


и т (К, ф,а }ехр[ва 24а, п (К,ф,а)= 


- 5 У (вал, (а, ) (3) 


п=ю ]= 
&2 = В 
+67 (К,а,@)В (а, А, (а, )ехр(тф ). 
Вид контура интегрирования Г в комплексной плоскости сх определяется использовани- 
ем принципа предельного поглощения при удовлетворении условиям излучения на бесконечности 


[9] и имеет следующий вид: обходит положительные особенности подынтегральной функции 
снизу, отрицательные - сверху, а на остальной части совпадает с вещественной осью. Функ- 


ции А,(@,@), В,(“,@) определяются при удовлетворении выражений (3) граничным усло- 


виям (2), С“ (В, а, ©}, С“ (В, а, 6), А, (©, @ ) и имеют довольно громоздкую структуру, 


ства нь) ата", 1 
оо)" (6 -а?)- съ(В,а, то 


например, 


С (вао)-"авута" сб®(В,а,в)=-а 6 (В, а, 6}; 
с85(В,а,в)= а? ‚(а с8°(В,а,0)=0; (4) 
Сова) = Ца" 

о, о итп. 


«Обезразмеривание» линейных координат проводится делением на величину а - внут- 


ренний радиус цилиндра. Соотношения, определяющие представления для С“? (К, ,@) полу- 
чаем из выражений для С“ (В, а ‚© ), заменой во всех выражениях функций Бесселя .Л,(7) на 


функции Неймана №, (7). 
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При подстановке представлений (3) в граничные условия (2) с использованием соотноше- 
ний закона Гука в цилиндрической системе координат, в преобразованиях Фурье по координате 
2 с разложением всех заданных и неизвестных функций в ряд Фурье по координате Ф получаем 
для каждого значения параметра И систему шести функциональных уравнений для определения 
неизвестных функций А; (а, ), В, («,@). Численные значения этих коэффициентов при за- 
данных значениях параметров проще получить, реализуя численный алгоритм их определения 
непосредственно на ЭВМ. Точность получаемых при этом результатов не ниже, чем при расчете 
по аналитическим представлениям. 

Итоговое НИ решений этой задачи имеет вид (3), где 


я, ( (в ‚ф,а > Ут Ф%® ( Ка, @ Хо (а )+ ФИ" (Ва, )Х (а Лехр(иф). (5) 


п=—ю = 


л + 
у(=) ь 1 ; 
Здесь Х’, @)=5-| [х/(е,ф)ехр[- Кпф+а 2 аф. 
—л —0 
При исследовании использована схема аналитико-численного алгоритма построения по- 
дынтегральных функций. В ней использованы аналитические представления (3), (4). Дальнейшие 
расчеты, связанные с решением систем линейных алгебраических уравнений для определения 


значений функций ФС (К,а,@) при заданных значениях параметров проводятся непосредст- 
венно на ЭВМ. 

Рассмотрим задачу об установившихся колебаниях бесконечного слоя с плоскопараллель- 
ными = границами т и - (- оо, оо}; -Н<2<0 при следующих граничных 
условиях: 

2=-Н: в,=У(х, уе", т, =У (уе, т, = (уе 


=0: о ‚= у (х,у)е“, = я (х,у)е*", ве у (х,у)е“. 


На бесконечности, как и ранее, задаем условия излучения энергии колебаний. Решение 
задачи строим методом интегральных преобразований с использованием принципа предельного 
поглощения [9]. В результате получаем следующее представление для амплитудной функции 


вектора смещения и, (х, у, 2) = {и\› и.» и, |: 


иуз(т,у,2) = ый [[22(.В,=)}ехрКу х+ В УмаВ- (7) 


(6) 


6 р 
Здесь У а 2)7,(у,В), 7=12,3; 


К=1 


бен +5 2)2(ейо 2-сйо (Н — 2 дей в „Н)+ 


и? А (и) 
+809 5йо „Н5йо (Н — =) 4и? (и? +0 (спо ,=-сйо ‚(Н -=ейоН)+ 
+(с в ,} (м? +в _ „Нэйо ›(Н-=)}; 


Г.(у,В,2)= (м? +а 2)? 2(спо 2-спо’ (Н -=ейо ,Н)+ 
ни а (и) 








+8и^о` ‚в ‚эй о „НУй д (Н-2)+ ди? (и? +0 2)? (спо 52 — СЙ О ‚(Н-=й @Н)+ 
+(о д ,) (и? +0 2 5й о Нэйо (Н-2) „... 





и? =у7+ В ?; с’=и?-0°; и; НЕЕ 2-Й; 
и р р 
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Д(и)=8(м? +02) (спо ‚Нео „Н —1)- 
5(? + 5] + 16и? (с та $ о по )5ло ,НзйоН 
РВ.) 71=1,2,3, получаем из 2. (7852), меняя во всех формулах 2 и Н-2 


местами, 


00 00 


ии, в) = [| [бе у)ехр(- {у х+ В ув. 


Полученные соотношения позволяют проводить расчет амплитудных характеристик век- 
тора смещения в слое при заданных на его гранях напряжениях У, (х,у )ехр(- 101). 


Рассмотрим исходную краевую задачу о возбуждении колебаний в конечном цилиндриче- 
ском слое в пространственной постановке (граничные условия (1)). В соответствии с принципом 
суперпозиции ищем в виде суммы решений вспомогательных краевых задач для бесконечного 
цилиндрического слоя с граничными условиями (2) и плоскопараллельного слоя (граничные усло- 
вия (6)): 

и(х, у,2) =и (х, у,2)+и.(К,ф,2). (8) 

Представления для компонент векторов и ‚ определены соотношениями (3),(4),(5) — для 


и и (7) для и,. Связь цилиндрических и декартовых координат определена соотношениями 
В? = +у’; ф=асю 2=2. При этом необходимо учитывать, что вектор смещений 
х 


И(х, у,2) =и(К,ф,2) может быть представлен в декартовых и цилиндрических координатах. 


Функции напряжений, определенные граничными условиями вспомогательных краевых 
задач, вообще говоря, неизвестны. Для их определения используем граничные условия исходной 
краевой задачи (1) и соотношения закона Гука в цилиндрических и декартовых координатах [6]. В 
результате получаем систему двенадцати интегральных уравнений относительно двенадцати не- 
известных функций напряжений х®(,2), У (х,у), (1=12,3; п=1,2,...,6) ‚ определяющих 


граничные условия ал задач следующей структуры: 














о (у, В, =) (у, В )ехру х" + Ву’) ав + 
— ХО (а ]ехр[а 2На= /®(2} )=12,3; т=1,2,3,...; 
ли Г /= 
Я АИ ехр[Кух" + Ву’) мв+ 


(9) 
о. И Дора ма = Ле) 





[=3 


Где Ф` («,,ф’, о )Х (а )+ 
НР 


-Ф( (“ ‚",ф”,@ ужо (а Лехр [2 2|4а+ 





а 
У (ху) + 
„(х,у) — 





| у (у, 8,0, 2) (у,В )ехр[Кух" + Ву’) |уаВ =Ф, (х,); 


41? и Го д 
п=1,2,3,4,5,6. 
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Здесь 


п 


1 (а) ных вет ) ао; 


пл 


р 1 + ; 
о а ео 


+00 


ХИ (е,ф) = У, Хо е)ехр[Итф }]; 


т=—<® 


х=Вс05ф; У = Ато; = чу; ф’=агю >. 
> 


Используя полученное решение исходной краевой задачи (граничные условия (1)), сфор- 
мулируем систему интегральных уравнений для реальной пространственной динамической кон- 
тактной задачи для поршневого пальца, находящегося в контакте с головкой шатуна. При поста- 
новке контактной задачи предполагается, что на части поверхности цилиндра заданы компоненты 
вектора смещения, а на остальной части — компоненты вектора напряжения. Рассмотрим случай, 


когда на внешней поверхности цилиндра А =а в области О: фсе [-ф о›Ф 92 [са [и,,^, | за- 
даны компоненты вектора смещения, например, и,(ф,2) =(И(2,ф). При контакте без трения ка- 
сательные напряжения в области <> полагаем равными нулю ( А = ТАМ = 0). Нормальные уси- 
лия в области контакта ( (В.Ф )) неизвестны и подлежат определению. Следовательно, в сис- 
теме (9) имеем 12 уравнений с 13 неизвестными. Для определения дополнительного неизвестного 
(нормального контактного давления (В.Ф )) используем условие равенства компоненты сме- 


щения поверхности цилиндра в области <> - и,(ф,2) смещению подошвы штампа ((2,ф). Это 
дополнительное уравнение имеет следующую структуру: 





2 : 2 ы 
ив = и (а,ф )- м2 (асозф,азтф, =)+и>,(а созф,азшф, 2) =0(2,Ф), (11) 
где и, определяется соотношениями (3) - (5); и›,и,> - соотношениями (7). В результате пра- 
вая часть (11) является интегральным оператором от неизвестных функций напряжений Х ьй,: 


Получаемое интегральное уравнение замыкает систему (9). Алгоритм получения представлений 
интегральных операторов достаточно прост при матричной записи. Однако при его реализации 
аналитические представления для подынтегральных функций операторов системы чрезвычайно 
громоздки и труднообозримы. Для исследования и построения решения систем предложено поль- 
зоваться численными методами, позволяющими производить расчет значений подынтегральных 
функций при фиксированных значениях параметров. Для заданных функций напряжений разра- 
ботаны алгоритмы и программы расчета интегральных представлений, определяющих операторы 
системы. 

Для разработки методики построения решения полученных выше систем интегральных 
уравнений необходим анализ свойств операторов систем. Исследование алгоритма построения и 
общих свойств операторов системы (9) подтверждает их непрерывность вдоль всей границы об- 
ласти, за исключением окрестности угловых точек. При этом можно показать, что в случае, когда 
нагрузка приложена на некотором удалении от этих точек, коэффициент при особенности реше- 
ния равен нулю [9]. При постановке контактной задачи дополнительное уравнение для определе- 
ния закона распределения контактных напряжений (11) является уравнением первого рода, свой- 
ства оператора которого аналогичны свойствам оператора соответствующей контактной задачи 
для бесконечного цилиндра. Следовательно, контактные напряжения имеют корневую особен- 
ность вдоль границы зоны контакта. 
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Анализ возможностей различных методов построения решения систем интегральных урав- 
нений при данных исходных определяет целесообразность использования комплекса методов по- 
строения ее решения для различных частотных диапазонов. В частности, для относительно низ- 
кочастотных колебаний (длина волны много больше размеров конечного цилиндра) можно доста- 
точно эффективно использовать аппроксимационный подход с использованием метода коллока- 
ций. При этом алгоритм построения приближенного решения систем сводится к следующему. Так 
как все неизвестные функции напряжений заданы в ограниченных областях, вводится их аппрок- 
симация в виде конечного отрезка разложения в ряд по ортогональной системе функций с неоп- 
ределенными коэффициентами и учетом порядка особенности (при ее наличии у данной функ- 
ции). При этом важно, чтобы ортогональная система функций была такова, что при подстановке 
аппроксимации в уравнение представлялась возможность точного вычисления наибольшего коли- 
чества интегралов в операторах, что повышает точность определения коэффициентов линейной 
алгебраической системы, получаемой при использовании метода коллокаций. 

Анализ особенностей строения данной системы определяет целесообразность использова- 


ния разложения функций Х», У, в степенной ряд по координатам К, 2 и в ряд по тригономет- 
рическим функциям (или ехр(тф)) по угловой координате. При определении функций 2,(2,ф) 


целесообразно выбрать аппроксимирующий полином следующей структуры: 


2(2,ф,0) =>, > С) ТФ %Ф)/ (#2 — 2Ф2—Ф?); Е=Ь2, 





где С,(@) - неопределенные постоянные; 7Т;(2), Т,(ф) - полиномы Чебышева. 


Подобное представление точным образом учитывает особенность контактных напряжений 
вдоль кромки штампа. Использование в разложении полиномов Чебышева определено их ортого- 


нальностью на требуемом отрезке с весом \/(#2 —=^) или /(ф?—ф?). 


Численные результаты и их обсуждение. Рассмотрим приложение разработанной выше чис- 
ленно-аналитической методики к уточненному расчету полей контактных напряжений в деталях 
цилиндропоршневой группы ДВС. Анализ проведем двумя способами: 1) при помощи метода ко- 
нечных элементов (использовался комплекс прикладных программ АМ$Уб) с учетом того, что об- 
ласть контакта заранее неизвестна и определяется в процессе вычислений; 2) решением системы 
интегральных уравнений с выделением математической особенности механических характеристик 
в нерегулярных точках области контактирования (область контакта постоянна). В качестве основ- 
ной расчетной модели был выбран ПП 406.1004020-01 (внешний диаметр Р=22мм, внутренний 
диаметр 14мм, длина [=64мм, масса л7= 0,113кг, материал - сталь 12ХНЗА) для двигателей 
марок 3МЗ-4061.10, 3МЗ-4062.10 и 3МЗ-4063.10, шатун в сборе 24-1004045-02 с втулкой 
406-1004052-10, поршень ЗМЗ 523-1004014. 

При реализации первого способа воспользуемся методикой работы [10]. В связи с тем, что 
точное распределение контактных напряжений между поршневым пальцем и бобышкой и между 
головкой шатуна и поршневым пальцем неизвестны и подлежат определению, в зоне контакта 
моделировали сетку контактных элементов (общим количеством до 36000). Взаимодействие дета- 
лей должно удовлетворять условиям контактирования, в соответствии с которыми точки одного 
тела не должны проникать внутрь другого. 

Для обнаружения условий проникновения задаются контактные узлы, а на поверхности 
поршневого пальца — контактные элементы. Вместе они образуют контактную пару. Контактный 
узел, попадая на поверхность контактного элемента, скользит вдоль нее с коэффициентом трения и. 
При вычислении контактных сил использовали известный метод штрафных функций. 

Результаты решения динамической задачи представлены в виде полей напряжений, де- 
формаций и перемещений деталей поршневой группы. На рис.1,2 показаны нормальные напря- 
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жения о,,о’,, соответствующие пятой моде колебаний и действующие в глобальных направле- 


х 
ниях У, Х в принятой системе координат. Видно, что в зоне контакта бобышки поршня с порш- 


невым пальцем возникает концентрация напряжений, обусловленная влиянием граничных эф- 
фектов. Поля напряжений имеют характерный пик в окрестности особой точки и ниспадающий 
характер на удалении. Похожая картина возникает в зоне контакта поршневого пальца с головкой 
шатуна. 





рмх =.109=—03 
эзмы =-.-450Е-+09 
эзмх =.1808-+09 
—-45090Е-+09 
—.-3з80овЕ-+09 
—-.3108Е-09 
—.240Е-+09 
—.170Е-+09 
—-998Е+08 
—.298Е-+08 
-ЗО2Е+ОЗ 
. 110Е-+09 
.180Е+089 














Рис. 1. Поле нормальных напряжений с, 


На рис. 1,2 видно, что кроме сжимающих напряжений, обусловленных контактным взаи- 
модействием, возникают значительные растягивающие напряжения, обусловленные изгибом и 
приводящие к овализации пальца. Попеременное воздействие этих напряжений приводит к появ- 
лению усталостных трещин и разрушению. В настоящее время в инженерной практике принято 
использовать приближенные зависимости, предложенные в работе [11]. Предполагается, что на- 
грузка по длине поршневого пальца распределяется по кусочно-линейному закону. Ясно, что та- 
кая схема нагружения будет давать весьма приближенные результаты расчета НДС пальца, так 
как не учитывает концентрации напряжений в особых точках, характерные черты которой «при- 
сутствуют» на рис.1,2. 





гмх =.109Е-ОЗ 

эмм =-.450Е+ОЗ 
5МХ =-361Е-+ОЗ 

.450Е+ОЗ 
.396Е-+ОЯЗ 
.-342Е+ОЗ 
.288Е+ОЯЗ 
-234Е+ОЗ 
.180Е+ОЗ 
-126Е+ОЗ 
.719Е-+ОЗ 
.179Е+08 
-361Е-+О08 


ЕТ 











Рис. 2. Поле нормальных напряжений о’, 
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При расчете по предлагаемой в данной работе численно-асимптотической методике пред- 
полагаем, что поле напряжений имеет корневую особенность вдоль границы зоны контакта. Это 
позволяет выделить сингулярные члены в разложениях механических характеристик и улучшить 
сходимость функциональных рядов для контактных напряжений. 

Для примера детально рассмотрим интенсивность локальной концентрации напряжений, 
возникающей в зоне контакта поршневого пальца с головкой шатуна (см. рис. 1,2) при нестацио- 
нарном динамическом нагружении (1). Подсчитаем относительное значение внутренней энергии 


Е [6], накопленной в нерегулярной зоне контакта <>. возле кромки головки шатуна. Значение 
Е представляет собой процентное отношение энергии, накопленной в области ©.„, к энергии, 
накопленной во всей области контакта головки шатуна с поршневым пальцем. Если взять начало 
координат в плоскости кромки головки шатуна, то область ©„ представляет собой цилиндриче- 
ский слой и задается неравенствами (координата 2 измеряется в миллиметрах): 
0,: фс Нл,ло2< [Еы. 

В таблице рассчитаны значения внутренней энергии на шести приведенных частотах, по- 
лученные при помощи метода конечных элементов (Ё\) и при помощи решения определяющей 
системы интегральных уравнений ( Е,). Величины Ё\ и Е, можно принять за показатели интен- 


сивности концентрации напряжений в области ©. 


Относительная внутренняя энергия, 
накопленная в области концентрации напряжений 






































Номер частоты эначение ее Е Е> 
колебаний, Гц 
1 30.7 51.3 56.4 
2 113.6 46.7 50.5 
3 326.2 46.1 47.7 
4 648.7 40.8 41.8 
5 864.5 38.8 38.5 
6 1023.6 38.0 34.9 





Из данных таблицы следуют несколько выводов: 
1. Для любой рассмотренной частоты колебаний наблюдается существенная концентра- 


ция напряжений в нерегулярной области контакта <>. , которая составляет лишь примерно 4% от 


общей площади контакта поршневого пальца с головкой шатуна. Это еще раз свидетельствует о 
невозможности применять в уточненных прочностных расчетах инженерные гипотезы [11]. 


2. Близость полученных значений Ё\ и ЕЁ, для каждого рассмотренного значения часто- 


ты говорит о том, что для оценки динамических характеристик контакта достаточно использовать 
предлагаемый в работе метод граничных интегральных уравнений и не привлекать сложный ко- 
нечно-элементный расчет, учитывающий изменение во времени области контакта. 

3. С увеличением порядкового номера моды колебаний интенсивность концентрации на- 
пряжений несколько уменьшается, приближаясь к постоянной величине. 

4. Результаты расчетов, проведенных по обеим методикам, практически совпадают, что 
свидетельствует о достоверности полученных данных и о возможном допущении постоянной об- 
ласти контакта в рассматриваемой области ©. 


5. Для достаточно больших значений частоты нагружения @« целесообразно применять 
конечно-элементные методики расчета, поскольку весь алгоритм решения определяющей систе- 
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мы интегральных уравнений контактной задачи базируется на том, что длина волны много боль- 
ше размеров поршневого пальца. 

Выводы. Анализ научных публикаций, посвященных построению аналитических подходов к ре- 
шению задач контактного динамического деформирования полых цилиндрических деталей маши- 
ностроения, показал, что в общем случае, когда изучение колебаний сводится к решению основ- 
ных граничных задач, возникают значительные математические трудности. Этим трудностям со- 
путствует существенное усложнение физической картины деформирования цилиндра. Кроме того, 
в области относительно высоких частот необходимо учитывать пространственный характер вол- 
новых движений частиц цилиндра. 

В связи с этим предложенная в работе методика обобщения метода суперпозиции для по- 
строения решения пространственной динамической задачи вибрационного деформирования по- 
лых цилиндрических тел обладает достаточным элементом научной новизны. Суть ее состоит в 
возможности представления общего решения в виде суммы решений вспомогательных задач для 
бесконечного цилиндрического слоя и бесконечного слоя с плоскопараллельными границами. 
Применение интегральных преобразований для решения вспомогательных задач позволяет полу- 
чить их решение в замкнутой форме. 

Большой интерес представляет оценка роли граничных условий в формировании спектра 
резонансных частот и форм колебаний полых цилиндрических деталей. С этой точки зрения в ра- 
боте рассмотрена динамическая смешанная контактная задача, которая имеет выраженное при- 
кладное значение. Получено интегральное уравнение, определяющее неизвестные напряжения в 
области контакта, которое замыкает определяющую систему интегральных уравнений относи- 
тельно неизвестных функций вспомогательной задачи. Таким образом, получены возможности 
численного исследования всех особенностей контактного деформирования полых цилиндрических 
деталей при общем случае нагружения их внешней поверхности. Рассмотрены возможности пре- 
небрежения в расчетах того факта, что область контакта ПП с головкой шатуна неизвестна и из- 
меняется в процессе деформирования. 

В качестве перспектив дальнейших исследований следует рекомендовать распространение 
предложенной численно-аналитической методики для анализа полей контактных напряжений для 
существенно анизотропных тел и тел с трещинами. 
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МУМЕКТСО-АМАТУТТСАЕ ТМУЕЗТТСАТТОМ ОЕ ОУМАМТС СОМТАСТ РВОСЕ$$Е$ 
ОЕ СУММОВТСАЕ МАСНТМЕ РАВТ$ 


Ма етайса| азрес5 оЁ {Пе зо!ибоп оЁ зраЧа! сопёасЕ упатс ргоБет оЁ сопасипа {йе сутанса! 
ЧефайЙ$ о сагз аге сопаегед. Урдайпта ог {Не те{о4 ог {Пе ирегро@оп 15 ргезетщеа, {Па а!о\м/5 © 
сопаисЕ апаУЧса! гезеагсй о{ {Не и\епзе-деГогтей сопа@оп апа {о таке ргасИса! гесоплтепааНоп$ 
ассог@та {о \еп$Ку оЁ оса! сопсепёгаНоп оЁ ргеззиге т итедуаг гопез оЁ сопёасЕ ПаЁ ргасйсаЙу 
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